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Claims] 

[Claim 1] A circuit simulation method comprising the steps of: 
5 (a) recognizing, from mask layout data for an integrated circuit, the shape of an 

electronic device to be analyzed which is provided in the integrated circuit, and obtaining 
data concerning the size of the electronic device to be analyzed; 

(b) determining the electrical characteristic of an electronic device for 
measurement, and measuring the size of each portion of the electronic device for 

10 measurement, as well as items each serving as an index of a stress applied to the electronic 
device to be analyzed; 

(c) extracting, based on at least the size of each portion of the electronic device for 
measurement, parameters from data concerning the electrical characteristic of the 
electronic device for measurement which has been determined in the step (b); and 

15 (d) utilizing a circuit simulator to select, from among the extracted parameters, a 

parameter suitable for each electronic device to be analyzed which is provided in the 
integrated circuit, and to perform circuit simulation in consideration of a stress applied to 
each electronic device to be analyzed. 

[Claim 2] The circuit simulation method of Claim 1, wherein in the step (b), at least an 
20 item serving as an index of a stress applied from an isolation insulating film to the 
electronic device to be analyzed is measured, and 

wherein in the step (d), the circuit simulation is performed in consideration of the 
stress applied from the isolation insulating film to the electronic device to be analyzed. 
[Claim 3] The circuit simulation method of Claim 1 or 2, wherein in the step (c), a 
25 plurality of parameters are extracted for each of the equal-sized electronic devices to be 
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analyzed, based on the items each serving as an index of a stress applied to the electronic 
device to be analyzed. 

[Claim 4] The circuit simulation method of any one of Claims 1 to 3, wherein the 
method further comprises, prior to the step (d), the step of inputting an additional model to 
5 the circuit simulator, the additional model being prepared based on measurement data that 
has been obtained in the step (b) and that serves as an index of a stress, and 

wherein in the step (d), a correction is made using the additional model when 
selecting a parameter suitable for each electronic device to be analyzed which is provided 
in the integrated circuit. 

10 [Claim 5] The circuit simulation method of any one of Claims 1 to 4, wherein the 
method further comprises, prior to the step (d), the step of preparing a reference table 
including pieces of information for associating each electronic device to be analyzed, 
which is provided in the integrated circuit, with the parameter that should be assigned to 
the electronic device to be analyzed, and the step of inputting the reference table to the 

15 circuit simulator, the reference table being prepared based on the items each serving as an 
index of a stress applied to the electronic device to be analyzed, and 

wherein in the step (d), the selection of the parameter suitable for each electronic 
device to be analyzed which is provided in the integrated circuit is automatically carried 
out using the reference table. 

20 [Claim 6] The circuit simulation method of Claim 5, wherein the reference table is 
used to associate each electronic device to be analyzed, which is provided in the integrated 
circuit, with a plurality of weighted parameters. 

[Claim 7] The circuit simulation method of any one of Claims 1 to 6, wherein the 
electronic device to be analyzed and the electronic device for measurement are each 
25 formed by a MIS transistor or a bipolar transistor. 
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[Claim 8] The circuit simulation method of Claim 7, wherein the electronic device to 
be analyzed and the electronic device for measurement are each formed by a MIS transistor 
comprising a gate electrode, a gate insulating film, an active region and an isolation 
insulating film surrounding the active region, and 
5 wherein the items, each serving as an index of a stress applied to the electronic 

device to be analyzed, include at least one of the position of the gate electrode in the active 
region, the size of the active region, and the width of the isolation insulating film. 
[Claim 9] The circuit simulation method of Claim 8, wherein the items, each serving 
as an index of a stress applied to the electronic device to be analyzed, further include at 
10 least one of the depth of the active region, a method for forming the isolation insulating 
film, the depth of the isolation insulating film, a material for use in forming the isolation 
insulating film, the size of the gate insulating film, and a material for use in forming the 
gate insulating film. 

[Claim 10] The circuit simulation method of Claim 8 or 9, wherein in the step (d), the 
15 circuit simulation is performed in consideration of a stress applied from the gate insulating 
film to the electronic device to be analyzed. 

[Claim 1 1] The circuit simulation method of any one of Claims 1 to 10, wherein in the 
step (b), at least an item that serves as an index of a stress applied from an interlayer 
dielectric film to the electronic device to be analyzed is measured, and 
20 wherein in the step (d), the circuit simulation is performed in consideration of the 

stress applied from the interlayer dielectric film to the electronic device to be analyzed. 
[Detailed Explanation of the Invention] 
[0001] 

[Technical Field to which the Invention Belongs] 
25 The present invention generally relates to circuit simulation methods to be 
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employed in designing a semiconductor integrated circuit. 

[0002] 
[Prior Art] 

Recently, in the field of LSIs such as semiconductor integrated circuits including 
5 MIS transistors, design specifications required for the integrated circuits have been more 
diversified and more complicated with increases in fineness of patterning for 
semiconductor elements, the number of semiconductor elements integrated on a chip, and 
the operating speed of each semiconductor element. 
[0003] 

10 In order to meet design specifications for various kinds of integrated circuits, each 

elementary circuit that has been designed or integrated circuit is subjected to circuit 
simulation so as to verify the function of each elementary circuit or the operation of the 
overall integrated circuit. In this case, parameters indicative of the characteristics of MIS 
transistors are extracted, and these parameters are used to predict how each MIS transistor 

15 operates. 
[0004] 

Normally, the obtainment of measurement data indicative of the characteristics of 
MIS transistors, utilized in the above-mentioned parameter extraction, requires the use of a 
semiconductor wafer on which several kinds of MIS transistors different in size (gate 
20 length L and gate width W) are formed. To be more specific, the principal characteristics 
of the MIS transistors formed on the wafer are measured, and parameters for the MIS 
transistors are extracted based on the electrical characteristics of the transistors. 
[0005] 

Hereinafter, parameters that have been used in conventional circuit simulation will 
25 be described in detail with reference to the accompanying drawings. 



4 



[0006] 

FIG. 12 is a graph showing the measurement results of drain current when different 
drain voltages (or source/drain voltages) Vd and gate voltages Vg are applied to a certain 
MIS transistor. From the measurement results shown in the graph, it can be understood 
5 that a drain current (Id)-drain voltage (Vd) curve is drawn for each of the gate voltages Vg 
(Vgl,Vg2and Vg3). 
[0007] 

In the conventional circuit simulation, the measurement values obtained with the 
drain current Id, drain voltage Vd and gate voltage Vg varied at appropriate steps are 
10 converted into Spice parameters, and these parameters are introduced into a circuit 
simulator. Furthermore, intermediate values between the measured points are interpolated 
using the Spice parameters, and are introduced into the simulator. 
[0008] 

FIG. 13 is a graph showing the relationship between the gate length L and drain 
15 current Id in the transistor when the drain voltage Vd and gate voltage Vg are kept 
constant. In the graph, "OD = 0.3|im" and "OD = 5.0^im" each represent the width of a 
source/drain region (active region) extending, at one side, from one end of a gate electrode 
to an isolation region in a direction parallel to the gate length. 
[0009] 

20 As can be seen from the characteristic curves obtained when "Id = Id 1" and "Id = 

Id 2" hold true, the characteristic of the transistor varies with a change in the gate length 
thereof. Therefore, measurement has to be also carried out with the transistor size (i.e., 
gate length L and gate width W) varied, and parameters responsive to the respective 
transistor sizes need to be provided based on the measurement. 

25 [0010] 



5 



However, it is actually difficult to provide a parameter for each transistor; 
therefore, a parameter is provided for several transistor sizes and is used in circuit 
simulation. 
[0011] 

5 FIG. 14 is a graph showing each range of transistor size to which a corresponding 

one of parameters divided into four is applicable. Specifically, illustrated in the graph is an 
example in which four parameters 1 to 4 are provided, and four transistor size ranges 1 to 4 
to which the corresponding parameters are applicable are provided. For example, circuit 
simulation is performed using the parameter 1 when the transistor size is in the range 1 in 
10 which a gate width is between Wl and W2 and a gate length is between L2 and L3, and 
circuit simulation is performed using the parameter 4 when the transistor size is in the 
range 4 in which a gate width is between W2 and W3 and a gate length is between LI and 
L2. 
[0012] 

15 FIG. 15 is a block diagram showing a conventional circuit simulation system. As 

shown in FIG. 15, a circuit simulator normally receives a netlist extracted from mask 
layout data, and parameters extracted from measurement values indicative of device 
characteristics. 
[0013] 

20 First, transistor size data 102, for example, is extracted from mask layout data 101 

including design information concerning a circuit to be analyzed, and the transistor size 
data 102 is inputted as a netlist 103 to the circuit simulator 100. As a matter of fact, the 
netlist 103 includes not only information concerning the transistor size but also 
information concerning capacitance and resistance. It should be noted that although FIG. 

25 15 shows the transistor data extracted from the mask layout data 101, data concerning 
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elements such as capacitor and resistor used to form a circuit is also actually extracted 
from the mask layout data 101. 
[0014] 

On the other hand, parameter extraction 105 necessary for circuit simulation is 
5 performed on measurement value data regarding a device for measurement (hereinafter, 
may also be called "device measurement data") 104, and the extracted parameters are 
inputted as parameters 106 to the circuit simulator 100. In the step of the parameter 
extraction 105, the obtained measurement value data 104 is converted into the parameters 
106. In the conventional method, not only the transistor size but also dopant concentration 
10 of a source/drain region and thickness of a gate insulating film, for example, have been 
considered. 
[0015] 

Next, the inputted parameters 106 are checked against the netlist 103 in the circuit 
simulator 100. Then, in the circuit simulator 100, an optimum model parameter 106a is 
15 selected for each transistor size 103a from among the inputted parameters 106, and circuit 
operation is simulated. 
[0016] 

For example, when a certain input signal is fed to the circuit to be analyzed, the 
simulation results indicating what kind of an output signal is obtained at an output terminal 
20 are provided as output results 107. In addition, circuit delay can be calculated in 
consideration of various resistances and capacitances. As the circuit simulator, a "SPICE" 
simulator or a tool obtained by making a modification to the simulator is generally used. 
[0017] 

Normally, circuit layout is modified with reference to the results of the simulation 
25 performed by the circuit simulator, and then simulation is performed again on the modified 
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layout by following the procedures similar to those described above. By repeating the 
procedures, an optimum circuit design can be carried out. 
[0018] 

[Problem that the Invention is to solve] 
5 In the above-described circuit simulation, based on design data regarding the 

transistor size and the inputted measurement data, the measurement data indicative of the 
electrical characteristic closest to the design size of each transistor is assigned to the 
corresponding transistor. Accordingly, it is basically impossible to eliminate an error 
between the calculated value obtained by the circuit simulation and the measurement value 
10 obtained using the actual circuit. Therefore, what is called for is to reduce the error 
between the calculated value in the circuit simulation and the measurement value to a level 
that causes no problem in circuit design. 
[0019] 

Suppose that the conventional method is performed using only a transistor size as a 
15 parameter when large design rules are used for an integrated circuit. Even in such a case, 
corrections are made in consideration of, for example, the shape of a gate electrode, the 
depth of a source/drain region and a dopant concentration thereof, thereby reducing an 
output error to a value that causes no problem from a practical standpoint. 
[0020] 

20 As the miniaturization of integrated circuits advances, however, the use of the 

conventional circuit simulation method has been causing an error between an actual circuit 
operation and an expected circuit operation to become more and more pronounced. Such 
an error in regard to circuit operation is aggravated when a MIS transistor or a bipolar 
transistor is provided, in particular, among various types of electronic devices. 

25 [0021] 
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It is expected that the miniaturization of integrated circuits continues to advance, 
and the use of design rules on the order of 0.13|im or less, in particular, strongly demands 
more precise and accurate circuit simulation. 
[0022] 

5 It is therefore an object of the present invention to provide more reliable and precise 

circuit simulation method applicable to integrated circuit design in which finer design rules 
are used. 
[0023] 

[Means of Solving the Problem] 

10 An inventive circuit simulation method includes the steps of: (a) recognizing, from 

mask layout data for an integrated circuit, the shape of an electronic device to be analyzed 
which is provided in the integrated circuit, and obtaining data concerning the size of the 
electronic device to be analyzed; (b) determining the electrical characteristic of an 
electronic device for measurement, and measuring the size of each portion of the electronic 

15 device for measurement, as well as items each serving as an index of a stress applied to the 
electronic device to be analyzed; (c) extracting, based on at least the size of each portion of 
the electronic device for measurement, parameters from data concerning the electrical 
characteristic of the electronic device for measurement which has been determined in the 
step (b); and (d) utilizing a circuit simulator to select, from among the extracted 

20 parameters, a parameter suitable for each electronic device to be analyzed which is 
provided in the integrated circuit, and to perform circuit simulation in consideration of a 
stress applied to each electronic device to be analyzed. 
[0024] 

According to the inventive method, the influence of stresses that has not been 
25 considered in conventional methods is factored into the parameters for the electronic 
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device to be analyzed, which have been provided for each size. Consequently, the circuit 
simulation can be performed accurately and precisely in consideration of a change in the 
characteristic of the electronic device caused by a stress applied thereto. 
[0025] 

5 In one embodiment of the inventive method, in the step (b), at least an item serving 

as an index of a stress applied from an isolation insulating film to the electronic device to 
be analyzed is preferably measured, and in the step (d), the circuit simulation is preferably 
performed in consideration of the stress applied from the isolation insulating film to the 
electronic device to be analyzed. In such an embodiment, all the stresses applied to the 
10 electronic device to be analyzed can each be approximated to the stress applied from the 
isolation insulating film. Accordingly, it becomes possible to relatively easily perform the 
accurate and precise circuit simulation in consideration of the stresses. 
[0026] 

In another embodiment of the inventive method, in the step (c), a plurality of 
15 parameters are preferably extracted for each of the equal-sized electronic devices to be 
analyzed, based on the items each serving as an index of a stress applied to the electronic 
device to be analyzed. In such an embodiment, the parameter that is close to the actual 
characteristic can be applied to each electronic device to be analyzed. As a result, it 
becomes possible to perform the circuit simulation with high degrees of precision, 
20 accuracy and reliability as never before. 
[0027] 

In still another embodiment of the inventive method, the method preferably further 
includes, prior to the step (d), the step of inputting an additional model to the circuit 
simulator, the additional model being prepared based on measurement data that has been 
25 obtained in the step (b) and that serves as an index of a stress. And in the step (d), a 
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correction is preferably made using the additional model when selecting a parameter 
suitable for each electronic device to be analyzed which is provided in the integrated 
circuit. In such an embodiment, even if the parameters extracted in the step (c) is not 
extracted in consideration of stresses, it is possible to perform the circuit simulation with a 
5 high precision in consideration of stresses. Besides, if the parameter extraction with 
stresses factored in is performed in the step (c), it is possible to further improve the 
preciseness and accuracy of the circuit simulation by using the additional model. 
[0028] 

In yet another embodiment of the inventive method, the method preferably further 
10 includes, prior to the step (d), the step of preparing a reference table including pieces of 
information for associating each electronic device to be analyzed, which is provided in the 
integrated circuit, with the parameter that should be assigned to the electronic device to be 
analyzed, and the step of inputting the reference table to the circuit simulator, the reference 
table being prepared based on the items each serving as an index of a stress applied to the 
15 electronic device to be analyzed. And in the step (d), the selection of the parameter 
suitable for each electronic device to be analyzed which is provided in the integrated 
circuit is preferably automatically carried out using the reference table. In such an 
embodiment, the time required for the simulation can be shortened. Therefore, such an 
embodiment is effective particularly when the number of the electronic devices to be 
20 analyzed is large. 
[0029] 

In one embodiment of the inventive method, the reference table is preferably used 
to associate each electronic device to be analyzed, which is provided in the integrated 
circuit, with a plurality of weighted parameters. In such an embodiment, since new 
25 parameters can be prepared by combining a plurality of parameters, the circuit simulation 
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can be performed with a higher degree of precision by using the new parameters. 
[0030] 

In another embodiment of the inventive method, the electronic device to be 
analyzed and the electronic device for measurement are each preferably formed by a MIS 
5 transistor or a bipolar transistor. Among various types of electronic devices, a MIS 
transistor and a bipolar transistor are likely to vary in electrical characteristic due to 
stresses applied thereto. Therefore, if parameters provided in consideration of stresses are 
used for a MIS transistor or a bipolar transistor, it becomes possible to easily perform the 
circuit simulation with a higher degree of precision as compared with the case where 
10 parameters provided in consideration of stresses are used for all types of electronic devices. 
[0031] 

In still another embodiment of the inventive method, the electronic device to be 
analyzed and the electronic device for measurement are each preferably formed by a MIS 
transistor including a gate electrode, a gate insulating film, an active region and an 

15 isolation insulating film surrounding the active region, and the items, each serving as an 
index of a stress applied to the electronic device to be analyzed, preferably include at least 
one of the position of the gate electrode in the active region, the size of the active region, 
and the width of the isolation insulating film. In such an embodiment, the influence of 
stresses can be factored in when the parameter extraction is performed, and furthermore, 

20 the influence of stresses can be factored in when the circuit simulation is performed. 
[0032] 

In yet another embodiment of the inventive method, the items, each serving as an 
index of a stress applied to the electronic device to be analyzed, preferably further include 
at least one of the depth of the active region, a method for forming the isolation insulating 
25 film, the depth of the isolation insulating film, a material for use in forming the isolation 
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insulating film, the size of the gate insulating film, and a material for use in forming the 
gate insulating film. In such an embodiment, the influence of stresses applied to the 
electronic device to be analyzed can be more closely reflected in the circuit simulation. As 
a result, the preciseness of the circuit simulation can be improved. 
5 [0033] 

In one embodiment of the inventive method, in the step (d), the circuit simulation is 
preferably performed in consideration of a stress applied from the gate insulating film to 
the electronic device to be analyzed. In such an embodiment, the influence of stresses 
applied to the electronic device to be analyzed can be more closely reflected in the circuit 
10 simulation. As a result, the preciseness of the circuit simulation can be improved. 
[0034] 

In another embodiment of the inventive method, in the step (b), at least an item that 
serves as an index of a stress applied from an interlayer dielectric film to the electronic 
device to be analyzed is preferably measured, and in the step (d), the circuit simulation is 

15 preferably performed in consideration of the stress applied from the interlayer dielectric 
film to the electronic device to be analyzed. In such an embodiment, again, the influence 
of stresses applied to the electronic device to be analyzed can be more closely reflected in 
the circuit simulation. As a result, the preciseness of the circuit simulation can be 
improved. 

20 [0035] 

[Mode for Carrying out the Invention] 

With an eye to improving the accuracy of circuit simulation, we conducted studies 
on factors that have not been considered in conventional circuit simulation, among factors 
that influence operations of electronic devices. After studying various kinds of factors, we 

25 found that the operation of a transistor is influenced by stresses applied from its 
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surroundings. 
[0036] 

Among stresses applied to a transistor, a stress applied from an isolation insulating 
film surrounding the transistor has the greatest influence on the transistor operation. The 
5 stress applied from the isolation insulating film formed by providing, for example, a 
shallow trench isolation (STI) region pressurizes or compresses an active region of the 
transistor. 
[0037] 

The characteristic curves "Id = Id 1" and "Id = Id 2" shown in FIG. 13 are the 
10 characteristic curves of MIS transistors that receive different stresses. The active regions 
of the transistors are different in size. Specifically, the characteristic curve "Id 1" is 
associated with "OD = 0.3|im", while the characteristic curve "Id 2" is associated with 
"OD = 5.0^im" ("OD" represents the width of a source/drain region extending, at one side, 
from one end of a gate electrode to an isolation region in a direction parallel to the gate 
15 length, and this width will be hereinafter called a "one-side OD width"). 
[0038] 

Suppose that the gate length is 0.3jim in FIG. 13. In that case, the drain current Id 
1 associated with "OD = 0.3nm" is about lSOpA/pm, and the drain current Id 2 associated 
with "OD = 5.0^m" is about 125^iA/|im. Accordingly, the drain currents Id 1 and Id 2 

20 differ from each other due to the different OD widths. From this fact, it can be seen that 
transistor characteristic is considerably influenced by a stress applied from an isolation 
insulating film. FIG 13 merely shows an example, and the electrical characteristic of a 
transistor varies depending on the conductivity type thereof, for example. However, it is 
true that the electrical characteristic of a transistor is considerably influenced by a stress 

25 applied thereto. 
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[0039] 

A stress applied from an isolation insulating film varies depending on the size of a 
transistor active region and/or a distance between the isolation insulating film and a gate 
electrode, for example. In light of this, the present inventors hit upon the idea of adding, as 
5 data to be measured, an active region size and/or a distance between a gate electrode and 
an isolation insulating film, for example, in order to utilize, as a new parameter for circuit 
simulation, an index of a stress applied to a transistor. 
[0040] 

Hereinafter, preferred embodiments of a circuit simulation method according to the 
10 present invention will be described with reference to the accompanying drawings. 
[0041] 
(First Embodiment) 

FIG. 1 is a block diagram illustrating a circuit simulation method according to a 
first embodiment of the present invention. According to the circuit simulation method of 
15 the present embodiment, a "SPICE" simulator or a modified one is utilized as in the 
conventional method, and an index of a stress applied to a transistor is used as a parameter 
in performing circuit simulation. 
[0042] 

As shown in FIG. 1, in the circuit simulation method of the present embodiment, 
20 netlist and parameter data are inputted to a circuit simulator. The netlist and the data are 
prepared as follows. 
[0043] 

The netlist 4 is extracted from mask layout data 1 regarding a circuit to be analyzed. 
[0044] 

25 To be more specific, the step of recognizing the shape of each transistor is 
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performed (hereinafter, called "transistor shape recognition" and identified by the reference 
numeral 2) based on the mask layout data 1. In the transistor shape recognition 2, each 
one-side OD width and each width of an isolation insulating film (hereinafter, may also be 
called an "isolation width") are recognized. 
5 [0045] 

Next, based on the results of the transistor shape recognition 2, the step of obtaining 
data including transistor size data 3a and transistor model recognition data 3b is performed. 
This step will be hereinafter called "data obtainment" and will be identified by the 
reference numeral 3. The transistor size data 3a obtained in this step includes pieces of 

10 information concerning transistor size (gate length and gate width), capacitance, resistance 
and wiring, for example. The transistor model recognition data 3b includes model names 
to be selected, which have been prepared manually based on each one-side OD width and 
each isolation width recognized in the transistor shape recognition 2. And the model 
names to be selected include data that serves as an index of stress. 

15 [0046] 

Then, the transistor size data 3a and transistor model recognition data 3b are 
inputted, as the netlist 4, to a circuit simulator 10. It should be noted that although not 
shown, not only data concerning transistor but also data concerning resistance and 
capacitance, for example, are actually inputted to the circuit simulator 10. 
20 [0047] 

Now, how the data for parameters 8 is prepared will be described. 

The data for parameters 8 is extracted from measurement values, i.e., device 
measurement data 5, which have been obtained using a device for measurement. The 
device for measurement is one selected or formed for the measurement, and is of the same 
25 type as the analyzed device. 
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[0048] 

Suppose that MIS transistors are used for the measurement. In that case, the size of 
each transistor is determined by its gate length L and active region width W, and the 
electrical characteristics of the MIS transistors different in size are determined, thus 
5 obtaining the device measurement data 5. Furthermore, the thickness of a gate insulating 
film, the shape of a source/drain region, a dopant concentration thereof, and a dopant 
concentration of a substrate, for example, are also measured under various conditions. In 
addition, in the present embodiment, factors related to stresses are also measured under 
various conditions. 
10 [0049] 

Subsequently, the step of recognizing the shape of each transistor (hereinafter, 
simply called "transistor shape recognition" and identified by the reference numeral 6) is 
performed based on the device measurement data 5. In the transistor shape recognition 6, 
measured one-side OD widths and isolation widths are recognized. 
15 [0050] 

Then, parameter extraction 7 is performed based on the results of the transistor 
shape recognition 6. Shown in FIG. 1 is an example in which parameter extractions 7a, 7b 
and 7c are carried out for three transistors that receive different stresses, based on 
parameters indicative of the stresses. Although the exemplary case where three kinds of 
20 stresses are applied is shown in FIG. 1, the parameter extraction may be carried out for four 
or more kinds of stresses. The parameter extraction 7 includes the step of converting the 
obtained device measurement data 5 into parameters 8 including model parameter groups 
8a, 8b and 8c provided in accordance with the magnitude of each stress. 
[0051] 

25 Next, the circuit simulator 10 receives the parameters 8 including the converted 



17 



model parameter groups 8a, 8b and 8c indicating the characteristics of the transistors each 
varied in accordance with the magnitude of stresses applied thereto. 
[0052] 

Then, upon receipt of the netlist 4 and the parameters 8 for the transistors, the 
5 circuit simulator 10 is utilized to select, based on the data included in the netlist 4, an 
optimum model parameter for each transistor size 4a from among the model parameter 
groups 8a, 8b and 8c provided in consideration of stresses, and circuit simulation is 
performed. In this case, information used for determining which of the model parameters 
is selected for each transistor is inputted based on the transistor model recognition data 3b. 
10 [0053] 

Subsequently, the parameters assigned to the respective transistors are used to 
obtain calculation results 9 from the circuit simulator 10. 
[0054] 

If the conventional circuit simulation method is performed in the circuit simulator 
15 10, a designer has no other choice but to assign identical parameters to equal-sized 
transistors that receive different stresses, since no parameters are provided in consideration 
of stresses in the conventional method. Therefore, an error is caused by characteristic 
variations resulting from different stresses, thus making it difficult to perform accurate 
circuit simulation. 
20 [0055] 

To the contrary, the circuit simulation method of the present embodiment makes it 
possible to select an optimum model parameter for each of equal-sized transistors, for 
example, from among the model parameter groups 8a, 8b and 8c in accordance with the 
stresses applied to the transistors. In the example shown in FIG. 1, an optimum model 
25 parameter can be selected for a transistor with a size of "Tr size 1" from among "Tr size la 
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model", "Tr size lb model" and "Tr size lc model" in accordance with the applied stress. 
[0056] 

Therefore, according to the circuit simulation method of the present embodiment, 
the preciseness and accuracy of the simulation are significantly improved as compared 
5 with the conventional method, and the simulation results can be utilized for circuit design 
in which finer design rules are used. Besides, according to the present embodiment, the 
number of stress-related factors to be measured, and the number of the parameter 
extractions are increased, thus making it possible to further improve the preciseness of the 
simulation. As described above, the circuit simulation method of the present embodiment 

10 is sufficiently adaptable to integrated circuit design in which finer design rules are used. 
Accordingly, the circuit simulation method of the present embodiment is preferably applied 
to circuit design in which design rules on the order of 0.13(im or less are used, for example. 
Naturally, the circuit simulation method of the present embodiment may be effectively 
used in designing already-existing integrated circuits. Consequently, with the use of the 

15 inventive circuit simulation method, innovative integrated circuits can be developed in a 
short period of time, and thus the products that meet the needs of the market can be 
provided without delay. 
[0057] 

The present inventors found items that should be measured in order to provide 
20 parameters in consideration of stresses, and these items are described below. 
[0058] 

Stresses applied to a MIS transistor include a stress applied from an isolation 
insulating film, a stress applied from a gate insulating film, and a stress applied from an 
interlayer dielectric film, for example. Among them, the largest one is the stress applied 
25 from an isolation insulating film. Therefore, at least the following items are each used as 
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an index for predicting the magnitude of the stress. 

• the size of an active region (length by width) 

• the length of the active region sandwiched between a gate electrode and an 
5 isolation insulating film (i.e., the position of the gate electrode in the active region) 

• the width of the isolation insulating film surrounding a transistor 

Hereinafter, the exemplary items to be measured will be specifically described with 
reference to the drawings. 
10 [0059] 

FIGS. 6 A and 6B are plan views each showing an exemplary MIS transistor 
including an active region and a gate electrode positioned in the active region. The 
transistors shown in FIGS. 6 A and 6B are equal in size, whereas their gate electrodes are 
different in position. Although not shown, each active region 61 is surrounded by an 
15 isolation insulating film (the same goes for FIG. 7). 
[0060] 

As shown in each of the plan views, a gate electrode 62 and dummy gate electrodes 
63 may be provided on one and the same active region 61 for manufacturing reasons, for 
example. In such a case, even if the transistors are equal in size, their electrical 
20 characteristics are different. It should be noted that the size of each transistor is 
determined by its gate length LI and active region width Wl. 
[0061] 

The electrical characteristic of each exemplary transistor varies depending on the 
position of the gate electrode 62 because a distance between the gate electrode 62 and the 
25 isolation insulating film is varied depending on the position of the gate electrode 62. In the 
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transistor shown in FIG 6A, the gate electrode 62 is located in an approximate center of 
the active region 61; on the other hand, in the transistor shown in FIG. 6B, the gate 
electrode 62 is located at one side of the active region 61 and adjacent to the isolation 
insulating film. Therefore, the gate electrode 62 of the transistor shown in FIG 6B is more 
5 susceptible to a stress applied from the isolation insulating film than the gate electrode 62 
of the transistor shown in FIG 6A, resulting in the transistors exhibiting different electrical 
characteristics. 
[0062] 

FIGS. 7 A through 7C are plan views each showing an exemplary MIS transistor 
10 including an active region and a gate electrode positioned in the active region. The active 
regions of the transistors are different in size, or the gate electrodes of the transistors are 
different in position. Shown in FIGS. 7 A through 7C are exemplary MIS transistors each 
having a gate length LI of 0.3 p.m and an active region width Wl of 10(im. Herein, "active 
region width" means the width of an active region extending in a direction parallel to the 
15 gate width. Furthermore, "active region length (one-side OD width)" herein means the 
width of an active region extending, at one side, from one end of a gate electrode to an 
isolation insulating film in a direction parallel to the gate length. 
[0063] 

FIG 7A shows the exemplary MIS transistor in which a gate electrode 60 is located 
20 in the center of an active region 64, and parts of the active region 64 located on both sides 
of the gate electrode 60 each have a length of 0.3 |im. 
[0064] 

FIG 7B shows the exemplary MIS transistor in which a gate electrode 60 is located 
in the center of an active region 65, and parts of the active region 65 located on both sides 
25 of the gate electrode 60 each have a length of 5.0|xm. 
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[0065] 

And FIG 7C shows the exemplary MIS transistor in which a gate electrode 60 is 
located at a left-side portion of an active region 66, and a part of the active region 66 
located on the left of the gate electrode 60 has a length of 0.3jxm while another part of the 
5 active region 66 located on the right of the gate electrode 60 has a length of lO.Ojxm. 
[0066] 

Since the MIS transistors shown in FIGS. 7 A and 7B are different in active region 
length, they receive different stresses from the isolation insulating films, and thus these 
MIS transistors exhibit different electrical characteristics. From this fact, it can be seen 
10 that the size of the active region can be used as an index of stress. 
[0067] 

Further, the entire width of the active region in the MIS transistor shown in FIG. 7B, 
extending in a direction parallel to the gate length, is almost equal to that of the active 
region in the MIS transistor shown in FIG 7C, extending in a direction parallel to the gate 
15 length; however, the gate electrodes of the MIS transistors shown in FIGS. 7B and 7C are 
different in position. Therefore, the gate electrodes of these transistors receive different 
stresses applied from the isolation insulating films, resulting in the transistors exhibiting 
different electrical characteristics. 
[0068] 

20 In view of the above, it is clear that the length of a part of an active region located 

on right of a gate electrode, and the length of a part of an active region located on left of a 
gate electrode can each be used as an index of stress. 
[0069] 

For example, in order to take into account the magnitude of each stress applied to 
25 the exemplary transistors shown in FIGS. 7 A through 7C, the parameter extractions 7a, 7b 
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and 7c are performed in accordance with the magnitude of each stress as shown in FIG 1 
in the present embodiment. Then, the parameters 8 including the results of the extractions, 
i.e., the model parameter groups 8a, 8b and 8c, are inputted to the circuit simulator 10, 
thereby making it possible to carry out the circuit simulation in consideration of the 
5 stresses. 
[0070] 

FIGS. 8 A through 8C are plan views each showing an exemplary MIS transistor 
including an active region 67 and a gate electrode 68. The transistors are surrounded by 
different-sized isolation insulating films. It should be noted that not only the active regions 

10 67 but also the gate electrodes 68 are similar in size and shape. To be more specific, the 
gate length of each gate electrode 68 is 0.3pm, the width of each active region 67 
extending in a direction parallel to the gate width is 10p,m, and the width of each active 
region 67 extending in a direction parallel to the gate length is 0.9pm (0.3pm + 0.3pm + 
0.3 jam). It should also be noted that the active regions 67 are equal in length and the gate 

15 electrodes 68 on the active regions 67 are similar in position. 
[0071] 

In the MIS transistor shown in FIG. 8A, an isolation insulating film 69 is formed to 
surround the periphery of the active region 67, and a semiconductor region (outward active 
region) 72 is formed to surround the periphery of the isolation insulating film 69. As 

20 shown in FIG. 8A, right and left portions of the isolation insulating film 69 located on the 
right and left of the active region 67, respectively, each have an isolation width of 4.0pm in 
a direction parallel to the gate length, while upper and lower portions of the isolation 
insulating film 69 located over and under the active region 67, respectively, each have an 
isolation width of 1 .0pm in a direction parallel to the gate width. 

25 [0072] 
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In the MIS transistor shown in FIG 8B, an isolation insulating film 70 is formed to 
surround the periphery of the active region 67, and a semiconductor region (outward active 
region) 73 is formed to surround the periphery of the isolation insulating film 70. As 
shown in FIG 8B, right and left portions of the isolation insulating film 70 located on the 
5 right and left of the active region 67, respectively, each have an isolation width of 4.0|im in 
a direction parallel to the gate length, while upper and lower portions of the isolation 
insulating film 70 located over and under the active region 67, respectively, each have an 
isolation width of 0.3|xm in a direction parallel to the gate width. 
[0073] 

10 In the MIS transistor shown in FIG. 8C, an isolation insulating film 71 is formed to 

surround the periphery of the active region 67, and a semiconductor region (outward active 
region) 74 is formed to surround the periphery of the isolation insulating film 71. As 
shown in FIG 8C, right and left portions of the isolation insulating film 71 located on the 
right and left of the active region 67, respectively, each have an isolation width of 0.3|im in 

15 a direction parallel to the gate length, while upper and lower portions of the isolation 
insulating film 71 located over and under the active region 67, respectively, each have an 
isolation width of 1 .Ojam in a direction parallel to the gate width. 
[0074] 

The right and left portions of the isolation insulating film 69 shown in FIG 8A and 
20 those of the isolation insulating film 70 shown in FIG 8B have identical isolation widths 
(i.e., 4.0jxm) in a direction parallel to the gate length. However, the upper and lower 
portions of the isolation insulating film 69 shown in FIG 8A and those of the isolation 
insulating film 70 shown in FIG 8B have different isolation widths (i.e., 1 .0|im in FIG 8A 5 
and 0.3 jim in FIG 8B) in a direction parallel to the gate width. In this case, the two 
25 transistors exhibit different electrical characteristics. This is because the magnitude of a 
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stress applied to a transistor varies in accordance with the isolation width of an isolation 
insulating film surrounding the transistor. 
[0075] 

Furthermore, the upper and lower portions of the isolation insulating film 69 shown 
5 in FIG 8A and those of the isolation insulating film 71 shown in FIG 8C have identical 
isolation widths (i.e., l.Ofim) in a direction parallel to the gate width. However, the right 
and left portions of the isolation insulating film 69 shown in FIG 8A and those of the 
isolation insulating film 71 shown in FIG 8C have different isolation widths (i.e., 4.0^m in 
FIG. 8A, and 0.3jxm in FIG. 8C) in a direction parallel to the gate length. In this case again, 
10 the two transistors exhibit different electrical characteristics. 
[0076] 

In view of the above, it is apparent that the size (isolation width) of an isolation 
insulating film surrounding a MIS transistor can be used as an index of a stress. 
[0077] 

15 FIGS. 9 A through 9C are plan views each showing another exemplary MIS 

transistor including an active region 67 and a gate electrode 68. The exemplary MIS 
transistors shown in FIGS. 9 A, 9B and 9C are surrounded by different-sized isolation 
insulating films 69a, 70a and 71a, respectively. The active region 67 and gate electrode 68 
of each MIS transistor shown in FIGS. 9A through 9C are similar to those of each MIS 

20 transistor shown in FIGS. 8A through 8C. And the isolation insulating films 69a, 70a and 
71a shown in FIGS. 9A through 9C, and the isolation insulating films 69, 70 and 71 shown 
in FIGS. 8 A through 8C have similar isolation widths not only in a direction parallel to the 
gate length but also in a direction parallel to the gate width. However, the MIS transistors 
shown in FIGS. 9 A through 9C are different from those shown in FIGS. 8 A through 8C in 

25 that semiconductor regions 72a, 73a and 74a located outwardly of the isolation insulating 
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films 69a, 70a and 71a, respectively, are each divided into four sections. In this case again, 
the magnitude of stresses applied to the MIS transistors shown in FIGS. 9 A through 9C are 
different from each other. 
[0078] 

5 In light of the above, the items, each used as an index of a stress indicative 

parameter, are summarized as follows. 
[0079] 

FIG 10 is a plan view of a MIS transistor, which is used to show exemplary main 
items that should be measured in order to obtain parameters in which the influence of 
10 stress is factored. Illustrated in FIG. 10 are an active region 75, a gate electrode 76, an 
isolation insulating film 77 and a semiconductor region (outward active region) 78. 
[0080] 

As shown in FIG 10, the main items each used as an index of stress in the circuit 
simulation of the present embodiment include, in addition to transistor size (gate length LI, 

15 and gate width Wl), the one-side OD widths of the inward active region 75, and the 
isolation widths of the isolation insulating film 77 surrounding the active region 75. To be 
more specific, the main items include: the one-side OD width (ODFL) of a left portion of 
the inward active region 75 located on the left of the gate electrode 76; the one-side OD 
width (ODFR) of a right portion of the inward active region 75 located on the right of the 

20 gate electrode 76; the isolation width (ODSL) of a left portion of the isolation insulating 
film 77 located on the left of the active region 75 in a direction parallel to the gate length; 
the isolation width (ODSR) of a right portion of the isolation insulating film 77 located on 
the right of the active region 75 in a direction parallel to the gate length; the isolation width 
(ODSU) of an upper portion of the isolation insulating film 77 located over the active 

25 region 75 in a direction parallel to the gate width; and the isolation width (ODSD) of a 
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lower portion of the isolation insulating film 77 located under the active region 75 in a 
direction parallel to the gate width. Herein, the widths ODFL and ODFR are collectively 
called "OD finger", while the widths ODSL, ODSR, ODSU and ODSD are collectively 
called "OD separate". 
5 [0081] 

FIGS. 11A and 11B are tables each showing a summary of items each used as an 
index of stress applied to the MIS transistor shown in FIG. 10. To be more specific, FIG. 
11B shows each index of stress applied to the MIS transistors shown in FIGS. 9A through 
9C. 
10 [0082] 

Each of the items shown above is measured as an index, and parameter extraction is 
performed based on the measurement, thus performing the high-precision circuit 
simulation using the extracted parameters in which the stresses applied to the MIS 
transistors are factored in. 
15 [0083] 

If the active region or isolation insulating film has a complicated shape, the other 
item having an influence on stress may be optionally added as the index. In such a case, 
the simulation can be performed with a higher degree of precision. 
[0084] 

20 Strictly speaking, a stress to be applied varies depending on the depths of isolation 

insulating film and active region, and a method for forming the isolation insulating film. 
Therefore, the circuit simulation can be performed with a higher degree of precision by 
taking into account the data concerning the depths of isolation insulating film and active 
region, and the method for forming the isolation insulating film. 

25 [0085] 
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Besides, a stress to be applied to a transistor varies depending on the material 
properties of an isolation insulating film. For example, there is a difference between a 
stress to be applied to a transistor in the case where SiC>2 containing no dopant is used, and 
a stress to be applied to a transistor in the case where BPSG (which is SiC>2 containing 
5 boron and phosphorus) is used. 
[0086] 

In addition, the size, thickness, material properties of a gate insulating film can be 
used as an additional index from the view point of stress. If an SOI substrate is used, the 
position of a buried oxide film, for example, can be used as an index of stress. 
10 Furthermore, by adding the thickness of an interlayer dielectric film as an index of stress, 
the simulation can be performed in consideration of the stress applied from the interlayer 
dielectric film. 
[0087] 

Although the circuit simulation method of the present embodiment has been 
15 described on the supposition that stress indicative parameters are assigned to MIS 
transistors, the parameters may also be assigned to bipolar transistors. In such a case, the 
items each used as an index of stress include: a distance between each of regions (which 
serve as a base, an emitter and a collector) and an isolation insulating film; and the size of 
the isolation insulating film. In addition, the present embodiment is also applicable to 
20 transistors other than those described above, capacitors, resistors, and diodes. The same 
goes for the embodiments described below. 

[0088] 
(Second Embodiment) 

FIG. 2 is a block diagram illustrating a circuit simulation method according to a 
25 second embodiment of the present invention. According to the circuit simulation method 
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of the present embodiment, an additional model is extracted from measurement data that 
serves as an index of the influence of stress, and is inputted to a circuit simulator. It should 
be noted that the same reference numerals as those used in the first embodiment are used in 
FIG 2 where appropriate. 
5 [0089] 

As shown in FIG 2, in the circuit simulation method of the present embodiment, 
not only a netlist 4 and parameters 8 but also an additional model 8d are inputted to a 
circuit simulator 10. The additional model 8d serves to correct the parameter assigned to 
each transistor in accordance with the magnitude of stress applied thereto. 
10 [0090] 

Parameter extraction 7A is performed on measurement values, i.e., device 
measurement data 5, each serving as an index of stress applied to a transistor, e.g., the OD 
finger, OD separate, and the depth of an isolation insulating film which have been 
described in the first embodiment, and the measurement values are converted into 
15 parameters and are inputted, as the additional model 8d, to the circuit simulator 10. 
[0091] 

Like the first embodiment, a netlist 4 is extracted from mask layout data 1 
regarding a circuit to be analyzed. That is, transistor shape recognition 2 is performed 
based on the mask layout data 1, and then the step of obtaining data including transistor 

20 size data 3a and transistor model recognition data 3b, i.e., data obtainment 3, is performed 
based on the results of the transistor shape recognition 2. The transistor size data 3a to be 
obtained in this step includes pieces of information concerning transistor size (gate length, 
and gate width), a dopant concentration of a source/drain region, capacitance, resistance 
and wiring, for example. The transistor model recognition data 3b includes model names 

25 to be selected, which have been prepared manually based on each one-side OD width and 



each isolation width recognized in the transistor shape recognition 2. And the model 
names to be selected include data that serves as an index of stress. 
[0092] 

In the method of the present embodiment, the step of recognizing the shape of each 
5 transistor, i.e., transistor shape recognition 6, is performed based on the size of each 
transistor as in the conventional method, and the parameter extraction 7 A is performed 
based on the measurement values, i.e., the device measurement data 5. Therefore, basically, 
one parameter is assigned to equal-sized transistors. 
[0093] 

10 However, in the circuit simulation method of the present embodiment, a correction 

is made using the additional model 8d in accordance with the magnitude of stress applied 
to each transistor when selecting a model parameter 8e for each transistor size 4a, thus 
making it possible to perform the simulation more precisely and accurately than the 
conventional simulation. In this embodiment, the selection of the parameter suitable for 

15 each transistor is carried out manually based on the prepared transistor model recognition 
data 3b. Alternatively, the selection may be carried out automatically utilizing computer 
software as in the embodiment described below. 
[0094] 

According to the method of the present embodiment, the additional model 8d that 
20 serves to correct parameters in accordance with the magnitude of stresses is added to the 
conventional model parameters 8e. Therefore, even if model parameters in which stresses 
are factored in are not available in a circuit simulator, the circuit simulation can be carried 
out with great precision in consideration of stresses by utilizing the additional model 8d, 
and thus high-precision output results 9 can be obtained. Furthermore, the preciseness of 
25 the simulation may be improved by preparing an additional model that indicates the 
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magnitude of stress in greater detail. 
[0095] 

In addition, the additional model can be also used when the parameter extraction is 
performed in consideration of the magnitude of stress as in the first embodiment. 
5 [0096] 

FIG. 3 is a block diagram illustrating a modified example of the circuit simulation 
method of the second embodiment. The modified example shown in FIG 3 is different 
from the second embodiment shown in FIG 2 in that parameter extractions are performed 
for equal-sized transistors in accordance with three magnitudes of stresses, for example. 

10 Furthermore, in a circuit simulator 10, model parameter groups 8f, 8g and 8h into which 
three additional models a, b and c are incorporated in accordance with respective applied 
stresses are prepared for equal-sized transistors. Therefore, an optimum model parameter 
can be selected from among the model parameter groups 8f, 8g and 8h for each of equal- 
sized transistors in accordance with applied stress. 

15 [0097] 

For example, although stress is factored into "Tr size la model" in the first 
embodiment shown in FIG. 1, stress is not factored into "Tr size la model" in the modified 
example of the second embodiment shown in FIG. 3. However, in the modified example of 
the second embodiment, the simulation can be performed in consideration of stress by 
20 making a correction using "additional model a". 
[0098] 

That is, in this modified example, the circuit simulation can be performed with a 
higher degree of precision since a correction is made to the three model parameter groups 
8f, 8g and 8h in consideration of stress by using the additional models a, b and c. 
25 However, for the additional models a, b and c, data that is more detailed than data used for 



the parameter extractions 7Ai, 7A2 and 7A3 has to be prepared. 
[0099] 

As described above, in the circuit simulation method according to the present 
embodiment or the modified example thereof, a correction is made in consideration of the 
5 influence of stress by using the additional model(s), thus making it possible to further 
improve the preciseness of the simulation. Consequently, the circuit simulation method 
according to the present embodiment can be sufficiently applied to circuit design in which 
finer design rules are used. 
[0100] 

10 (Third Embodiment) 

FIG 4 is a block diagram illustrating a circuit simulation method according to a 
third embodiment of the present invention. It should be noted that the same reference 
numerals as those used in the first embodiment are used in FIG. 4 where appropriate. 
[0101] 

15 The circuit simulation method of the third embodiment is different from that of the 

first embodiment in that the method of the third embodiment utilizes a reference table 12 
for associating each transistor size 4a with an optimum model parameter selected from 
among model parameter groups 8a, 8b and 8c. 
[0102] 

20 In the first embodiment, for the selection of a model parameter most suitable for 

each transistor size 4a included in the netlist 4, information used for associating each 
transistor size with each model parameter is manually inputted to the transistor model 
recognition data 3b by a designer. To the contrary, in the circuit simulation method of the 
third embodiment, the netlist 4, data for parameters 8, and reference table 12 are inputted 

25 to a circuit simulator 10. In this case, in the transistor model recognition data 3b, only one- 
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side OD widths and isolation widths are inputted, and no model names are inputted unlike 
the first embodiment. In the circuit simulator 10, a model parameter suitable for each 
transistor size 4a is automatically selected from among the model parameter groups 8a, 8b 
and 8c based on information provided in the reference table 12. 
5 [0103] 

After transistor shape recognition 2 performed using mask layout data 1 and 
transistor shape recognition 6 performed using device measurement data 5, a transistor 
reference table 11 is manually prepared based on the results of both of the transistor shape 
recognition 2 and the transistor shape recognition 6. And the prepared transistor reference 
10 table 11 is automatically inputted, as the reference table 12, to the circuit simulator 10. In 
the reference table 12, for example, "Tr 1" is associated with a parameter "Tr la", and "Tr 
2" is associated with a parameter "Tr 2b". 
[0104] 

According to the present embodiment, in the circuit simulator 10, the reference 
15 table 12 is utilized to automatically select a model parameter that is most suitable for each 
transistor size. Therefore, even if the number of transistors is increased, it does not take 
much time to analyze the transistors. This is because, although the time required for the 
preparation of the reference table 12 does not change much with an increase in the number 
of transistors, the time required for the analysis performed by the circuit simulator is 
20 shorter as compared with the case where the analysis is performed manually. 
[0105] 

Consequently, according to the circuit simulation method of the present 
embodiment, if the number of transistors is large, the time required for the analysis can be 
shorter as compared with the first embodiment. The preciseness of the simulation in the 
25 present embodiment is similar to that of the simulation in the first embodiment. 
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[0106] 

The present embodiment has been described as an example in which the reference 
table is utilized in the method of the first embodiment. Alternatively, the reference table 
may be effectively utilized when the additional model is used as described in the second 
5 embodiment. 
[0107] 
(Fourth Embodiment) 

FIG 5 is a block diagram illustrating a circuit simulation method according to a 
fourth embodiment of the present invention. It should be noted that the same reference 
10 numerals as those used in the third embodiment are used in FIG. 5 where appropriate. The 
fourth embodiment is different from the third embodiment in that a transistor reference 
table 13, a combined reference table 14, and a combined model parameter group 8A are 
added in the fourth embodiment. 
[0108] 

15 As shown in FIG 5, according to the circuit simulation method of the present 

embodiment, a plurality of parameters can be selected for one transistor by utilizing the 
combined reference table 14 in a circuit simulator 10. 
[0109] 

The circuit simulator 10 receives a netlist 4, model parameter groups 8a, 8b and 8c, 
20 and the combined reference table 14 that has been prepared in advance based on the 
transistor reference table 13. In the present embodiment, the combined reference table 14 
is utilized to select a plurality of model parameters for one transistor, and the circuit 
simulation is performed using the combined model parameter group 8A determined by the 
weighting of each model parameter. In this manner, output results are obtained. 
25 [0110] 



In the example shown in FIG 5, model parameters "Tr la" and "Tr lb" are selected 
for a transistor "Tr 1" by utilizing the combined reference table 14, and each of the 
parameters is weighted. For example, if a stress applied to the transistor "Tr 1" is at an 
intermediate value exactly between the value of the parameter "Tr la" and the value of the 
5 parameter "Tr lb", fl model "fl (Tr la, Tr lb) = (Tr la * 0.5 + Tr lb x 0.5)" is assigned to 
the transistor "Tr 1". Thus, if a stress applied to a transistor is at an intermediate value 
between the values of the model parameters included in the model parameter groups 8a, 8b 
and 8c obtained by the parameter extractions 7a, 7b and 7c, a combined model parameter 
indicative of the intermediate value can be prepared and assigned to the transistor. In the 

10 third embodiment, however, an optimum model parameter indicative of the value of a 
stress is selected only from among the model parameter groups 8a, 8b and 8c which have 
been obtained by the parameter extractions 7. To the contrary, according to the fourth 
embodiment, the circuit simulation can be performed using a combined model parameter 
indicative of an intermediate value between the values of model parameters, thus making it 

15 possible to obtain high-precision output results. 
[0111] 

As described above, according to the circuit simulation method of the present 
embodiment, a plurality of parameters are selected for one transistor utilizing the combined 
reference table 14, and a combine model parameter is newly generated based on these 

20 parameters. As a result, it becomes possible to further improve the preciseness and 
accuracy of the circuit simulation. Which of the parameters is selected for a certain 
transistor and how the extracted model parameters are weighted may be determined by 
taking into consideration each index of stress such as the shape of an active region and the 
position of a gate electrode. 

25 [0112] 
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In the circuit simulation method of the present embodiment, two parameters do not 
necessarily have to be selected for one transistor, but three or more parameters may be 
selected for one transistor. 
[0113] 

Furthermore, the circuit simulation method of the present embodiment may be 
effectively used when the additional model is utilized as in the second embodiment. 

[0114] 
[Effect of the Invention] 

According to the circuit simulation method of the present invention, the influence 
of stress applied to the electronic devices is factored in the parameters, thereby making it 
possible to improve the preciseness and accuracy of the circuit simulation. Thus, design of 
the integrated circuit which has been increasing in fineness is carried out swiftly and new 
products are introduced into the market in a short time. 
[Brief Explanation of the Drawings] 

[FIG. 1] 

FIG. 1 is a block diagram illustrating a circuit simulation method according to a 
first embodiment of the present invention. 
[FIG. 2] 

FIG. 2 is a block diagram illustrating a circuit simulation method according to a 
second embodiment of the present invention. 
[FIG. 3] 

FIG. 3 is a block diagram illustrating a modified example of the circuit simulation 
method according to the second embodiment. 
[FIG. 4] 

FIG. 4 is a block diagram illustrating a circuit simulation method according to a 
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third embodiment of the present invention. 
[FIG. 5] 

FIG. 5 is a block diagram illustrating a circuit simulation method according to a 
fourth embodiment of the present invention. 
[FIG. 6] 

FIGS. 6A and 6B are plan views each showing an exemplary MIS transistor 
including an active region and a gate electrode positioned in the active region. The 
transistors are equal in size, whereas their gate electrodes are different in position. 

[FIG. 7] 

FIGS. 7A through 7C are plan views each showing an exemplary MIS transistor 
including an active region and a gate electrode positioned in the active region. The active 
regions of the transistors are different in size, or the gate electrodes in the active regions 
are different in position. 

[FIG. 8] 

FIGS. 8 A through 8C are plan views each showing an exemplary MIS transistor. 
The transistors are surrounded by different-sized isolation insulating films. 
[FIG. 9] 

FIGS. 9A through 9C are plan views each showing another exemplary MIS 
transistor. The transistors are surrounded by different-sized isolation insulating films. 
[FIG. 10] 

FIG. 10 is a plan view of a MIS transistor, which is used to show exemplary main 
items that should be measured in order to obtain parameters in which the influence of 
stresses are factored in. 

[FIG 11] 

FIGS. 11A and 11B are tables each showing a summary of items each used as an 
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index of stress applied to the MIS transistor shown in FIG 10. 
[FIG. 12] 

FIG. 12 is a graph showing the electrical characteristics of a MIS transistor having 
a certain size when different gate voltages Vg are applied. 
5 [FIG. 13] 

FIG. 13 is a graph showing the relationship between gate length and drain current 
in the transistor when the drain voltage Vd and gate voltage Vg are kept constant. 
[FIG. 14] 

FIG. 14 is a graph showing exemplary transistor size ranges to each of which a 
10 corresponding one of parameters for circuit simulation is applied. 

[FIG. 15] 

FIG. 15 is a block diagram showing a conventional circuit simulation system. 
[Explanation of Reference Numerals] 
1 Mask layout data 
15 2,6 Transistor shape recognition 
3 Data obtainment 
3 a Transistor size data 

3 b Transistor recognition data 

4 Netlist 

20 4a Transistor size 

5 Device measurement data 
7, 7a, 7b, 7c Parameter extraction 

7 A, 7Ai, 7A 2 , 7A 3 Parameter extraction 
8 Parameter 
25 8a, 8b, 8c Model parameter group 
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8d Additional model 
8f, 8g, 8h Model parameter group 
8e Conventional model parameter 
8A Combined model parameter group 

9 Output results 

10 Circuit simulator 
11,13 Transistor reference table 
1 2 Reference table 

14 Combined reference table 

60, 62, 68 Gate electrode 

61, 64, 65, 66, 67 Active region 
63 Dummy gate electrode 

69, 69a, 70, 70a, 71,71a Isolation insulating film 
72, 72a, 73, 73a, 74, 74a Semiconductor region 

75 Active region 

76 Gate electrode 

77 Isolation insulating film 

78 Semiconductor region 
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[Name of the Document] ABSTRACT 
[Summary] 

[Purpose] A purpose of the invention is to provide a circuit simulation method used 
for designing integrated circuits under finer design rules with improved preciseness and 
5 accuracy. 

[Solution] In an inventive circuit simulation method, simulation is performed utilizing 
a circuit simulator, based on a netlist prepared using mask layout data for a circuit, and 
parameters obtained from measurement data concerning the characteristic of each 
transistor. The parameters are extracted from the measurement data based on not only the 
10 transistor size but also a stress applied to the transistor. Therefore, the circuit simulation 
can be performed with precision and accuracy never before possible, in consideration of a 
change in the characteristic of the transistor which is caused by the stress applied thereto. 
[Selected figure] FIG. 1 
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